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Abstract 

We show that an high temperature expansion at fixed order parameter can be derived for the 
quantum Ising model. The basic point is to consider a statistical generating functional associated 
to the local spin state. The probability at thermal equilibrium of this state refiects directly the 
occurrence of a spontaneous symmetry breaking. It is possible to recover the expansion around 
the mean field in the system dimensionality if the "direction" in the Hilbert space of local spin 
states is suitably chosen. Results for the free energy at the critical temperature, as a function of 
the transverse field, in first order approximation in the inverse system dimensionality are compared 
with those of the standard approach. 
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The recent interest in decoherence[0 phenomena in connection with quantum computing 
leads naturally to the study of a local spin state interacting with an environment. It has 
been shown that quantum computation can be performed acting on a nuclear spin at thermal 
equilibrium 10]. In this perspective it is interesting to study the simplest quantum spin 
model, i.e the Ising model in a transverse field. This model has an exact solution in one 
dimension]^], and approximate techniques, such as the random phase approximation |Q and 
exact cumulant series expansion § have also be applied for higher order. The main result 
is due to Suzuki's work0 who showed the correspondence with a classical model in d+1 
dimensions. A systematic expansion around the mean field result is however still lacking. 
As we shall see such an expansion can be introduced if the phase transition is characterized 
by a statistical generating functional associated to the local spin state. The probability of 
finding the spin in the i-th site in the general state \si >= cos Ti> +e*'^' sin yl ii> when 
the system is in thermal equilibrium is 



Tr 



Si >< Si 



It is convenient to express the projection operators in terms of the Pauli spin operators 

\si >< Si\ = ^ + Si- rii (2) 

Here rii (sin 6i cos v^i, sin Oi sin ipi, cos Oi) and Si (S'f , , S'f ) are respectively the unitary vector 
which defines the symmetry breaking direction, and the spin - | vector operator. Taking 
into account Eq. (0), the local spin state probability becomes 

(a \ ^ Tr [ex^{-(3H)Si ■ Ui] 

Note that the non trivial local spin probability Pi — \ corresponds to the magnetization in 
the symmetry breaking direction. It is immediately seen that pi is directly related to the 
symmetry breakdown. In the quantum Ising model 

H = -Y.J^.sts;-Y.h,s^ (4) 

the only symmetry is the spin inversion S^ ^S^. Above the critical temperature Tc, the 
magnetization is directed along the transverse field direction, while in the ordered phase 
a non vanishing component in the z direction appears. The local state probability of the 
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symmetrical phase (above Tc) is 

where n±i is the unitary vector orthogonal to the z direction. When the symmetry is broken 
(below Tc) 

^SB(n ,^ ^ _ ^S(n ,^ ^ , Tr [ex^{- [3 H)Sl COS Oj] 
P. (e.¥^.)-pa^.¥'.)+ Tr[eM-m] 

We see that the qualitative change of the local probability is directly related, as expected, to a 

non vanishing statistical average of S^. The non trivial part of the local spin state probability 

is the first moment of the statistical generating function G in the limit of vanishing source 

field A,-. 



G = \n{Tr 



■ rii 



(7) 



We see that for vanishing Aj, G becomes the free energy. Moreover in the absence of a 
transverse field and for 6i = the statistical functional G reduced to that of the Ising model 
in the presence of a symmetry breaking field along the z - direction hi = ^. 

Expansion around the mean field are usually derived by introducing a generalization of 
the free energy in the presence of symmetry breaking fields and then by Legendre transform 
which makes the order parameter the independent variable of the problem. This general- 
ization is usually accomplished introducing symmetry breaking fields as perturbation in the 
system Hamiltonian. We consider here an alternative generalization considering the gener- 
ating function of Eq. (0). The main advantage of the present approach is the possibility of 
a straightforward expansion in the inverse temperature. At infinite temperature there is a 
complete degeneracy of the local spin state, which is removed only for the presence of the 
symmetry breaking fields. At finite temperature, the correction due to spin - spin interaction 
and the transverse field compete each other to determine the local spin state probability and 
the local magnetization. 

Let us to introduce the Legendre transform of the generating function W{m) = G{X) — 
Here rrii =< Si ■ rii > and angular brackets stand for 



^ Tr 
<0> 



e-^^^exp{J:^HP){S^■n,)) 



(8) 



Tr [e-P^exp{E^KmSi■n,))] 

Note that rrii, as a consequence of the relation between rrii and the generating function 
rrii = become the independent variables of W. On the other hand, the source field Aj 



depends on rrii and (3 



dW 

drrii 



(9) 



The free energy in the hmit of vanishing source field is given by W calculated for the 
"magnetization" which satisfies extremum condition |^ = 0. This is the same procedure 
which applies to the free energy of the Ising mo del Following reference (2), we exploit the 
limit of vanishing /3 to fix the relation between and Aj(0) 
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The next step is to derive the high temperature expansion of W. The expansion up to the 
second order in j3 is related to the following quantities. 



W\f3=o 

dW 



d/3 



(5=0 



13=0 



^In 2cosh(^M2)^ 

>o 

< >n -(< H >n 



Xi{0)mi 



(11) 



s2 _r>sr ^ 



< H{Si ■ Hi - rrii) >o + 



13=0 



E 



dp 



13=0 



dp 



< {Si ■ rii - mi){Sj ■ Uj - rrij) >o 



(3=0 



There the symbol <>o stands for the statistical average defined in Eq. evaluated at 
P = 0, and the derivative of Aj with respect to P are obtained from Eq. @ 



dp 



(3=0 



dpdrrii 



(12) 



(3=0 



A straightforward calculation gives 



— m, 



+P ^ Jijuiirrij cos 9i cos 6'j + /? ^ himt sin 9i cos ipi 



+ — E('^*J + Jj,i){Ji,k + Jk,i)mjmk cos 6k cos 6j sin^ + — ^ h'^{l - sin^ 9i cos^ Lpi) 
— — E('^«J + COS 6'j sin 6'j cos 6'j cos ipi — — E('^*j + "^i.O^^^i "^o^^ 6'^ sin^ 6i 
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(13) 



4 



Eq. ([lOD has been used to eliminate the dependence from Aj(0). Note that the first fine 
of the left hand side of Eq. (|1^) corresponds to the mean field approximation. In the 
limit of vanishing transverse field, i.e. when the symmetry for rotation around the z axis is 
recovered, the dependence on ipi disappears as expected. Moreover if the symmetry breaking 
fields is fixed along the z - axis, we recover the free energy expansion of the Ising model[Q. 
This direction is however selected from the extremum condition on the free energy which 
corresponds to a minimum in mean field approximation. In such a case it is easily shown 
that, considering a constant interaction among next neighbors sites, which scales with the 
system dimensionality (Jjj + Jj^i = J/D), the high temperature expansion becomes an 
expansion in This result is actually a consequence of the compensation of contributions 
which are of second order in /? and zero order in As we shall see this compensation is 
valid also in the presence of transverse field with a suitable choice of the symmetry breaking 
fields direction. 

Let us consider for simplicity the case of an homogeneous symmetry breaking = m, 
9i = 9 and (pi = (p. The minimum of the free energy in ip gives, as expected = 0. The 
probabilistic functional becomes 



+ m ) In ( — \- ni 



W{m,9) = N 

P ^— cos^ 9 + hm sin 9^ + 



m I In 



+ 



— cos 9 



[Jmsm9 — h) — sm 9 



D 



+ 



m 



cos^ 9 + 



16 



cos 



(14) 



Considering only contributions in (3 up to the first order we recover the well known results 
of the mean field approximation, by an extremum condition in m and 9. The first condition 

dW n ■ 

= gives 

am c3 



m 



- tanh 
2 



— (^Jm cos^ 9 + h sin 9^ 



(15) 



corresponds to the vanishing of the symmetry breaking fields, i.e. to the limit in which the 
probabilistic generating function W reduces to the free energy. The second corresponds to 
the selection of the direction of the symmetry breaking field corresponding to an extremum 
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of free energy = 0. 

cos^ {JmsinO — h) = (16) 

It is important to note that the first order expansion in j3 plus the minimum condition in 
6, if gives the exact solution associated to the mean field Hamiltonian. It means that the 
rotation in the 6 and space that diagonalizes the mean field Hamiltonian coincides with 
the choice of 6 and (p which minimize the first order correction to the energy. 

As far as the second order terms in f3 are concerned, we note that the "spurious" term 
which does not scale with the inverse dimensionality D, can be neglected close to the mag- 
netization directions selected in the first order approximation (sin = under the critical 
temperature, and, cos^ = above). The solution for the extremum condition, below the 
critical temperature for 6 and m modifies as follows: 



m 



i tanh [^Jm-^ {^h^ + J\l - 4m^)) } (18) 

Within this approximation we find an expansion in ^ which is the generalization of the 
Ising model results[|]. In the limit of vanishing transverse field we obtain, as expected, the 
Ising model results. 

The critical temperature is identified as the temperature at which the symmetry breaking 
direction coincides with that of the transverse field, i.e sin 6' = 1. Solving for the extremum 
conditions we obtain a self consistent equation for the critical temperature. 

= '^arctanh {^—^^^^ | (19) 



The solution for Eq. (19) is defined in a limited range of dimensionality and transverse field 
amplitude. The lowest critical temperature and highest transverse field are given by 

2h 

= tanh 




= ^ l-2lK^ + 77^ (20) 



In order to compare our results with those derived in reference (9), we must take into account 
that this model is the same of Eq. (|D with spin operators substituted by Pauli operators. 
This implies the following scaling in the Eqs. (p^ (p!9D:m ^ ^, J —>■ 4 J, h —>■ 2h. 



-1 



-1.5 





D=10 




D=10 (Stratt) 




D=5 




D=5 (Stratt) 




\ \ 


\ 

\ 

\ 

\ 


\ \ 


\ 

\ 

1 

1 

1 
1 


'■\ \ \ 
\1 \ 

f. 



0.2 0.4 0.6 0.8 1 

Transverse Field h/J 



FIG. 1: The free energy for site at critical temperature as function of transverse field 
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FIG. 2: The critical temperature as function of transverse field. The lower curve gives the temper- 
ature and the corresponding transverse field which determines the limit of our approximation. 

In Fig. (1) the free energy at critical temperature obtained in the framework of the 
standard approach of reference (9) is compared with our results. We find an improvement 
for /i 7^ 0. The same comparison for the critical temperature as function of the transverse 
magnetic field is performed for dimensionality greater than four in Fig. (2). It is easily 
seen that the free energy is lowered because of corrections of order in the direction 9 
of the magnetization. The presence of this extra variational parameter seems to be the 
main advantage of the present approach. As a concluding remark we would like to comment 
about possible extensions of this approach to other quantum system. We think about system 
defined on a lattice where in the energy there is a "local" term, which depends only on the 
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degrees of freedom associated to a single site, plus an interaction which takes into account the 
coupling between degrees of freedom of different lattice sites. If for a particular choice of the 
system parameters there is a symmetry in the local energy operator we can treat interaction 
and symmetry breaking part of the local Hamiltonian as a perturbation. The degeneracy of 
the system in the limit of vanishing perturbation, can be removed by the introduction of the 
source field, which determines the "orientation" in the Hilbert space of the system ground 
state. In the case of the quantum Ising system the local energy is zero and then any local 
spin state is allowed. The source field A determines the local state. The expansion parameter 
is associated to the inverse temperature f3. Taking into account the perturbation we can 
develop an expansion whose first term is the mean field approximation. The convergence 
of this expansion can be improved if a direction selection mechanism, analogous to that 
discussed previously applies. It means that the minimum of the first order approximation 
of the free energy (mean field) determines a direction in the degenerate state Hilbert space 
and next order contributions give only small correction to that direction. We are referring 
to lattice models where the local energy is determined by the site particle occupation and 
particle - particle interaction on the same site, and the interaction among sites is associated 
to hopping of particles from on site to next neighbors site (Hubbard model). For particular 
values of the chemical potential there is a degeneracy of the statistical weight of the two 
lowest energy states of the local Hamiltonian. We use the subspace associated to these 
states to introduce a general local state as a superposition of the two independent ground 
states. The spontaneous symmetry breakdown will be discussed in terms of the probability 
of a local ground state at thermal equilibrium. Application to repulsive Hubbard model for 
bosons on a lattice, and attractive Hubbard model for fermions are in progress. Preliminary 



results concerning the boson system can be found in reference 1 10 
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